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NJL-type effective models represent a low-energy realization of QCD 
and incorporate pertinent aspects such as chiral symmetry and its sponta- 
neous breaking, the center symmetry in the heavy-quark limit as well as the 
axial anomaly. One such model, the Polyakov-quark-meson model for three 
light quark flavors, is introduced in order to study the phase structure of 
strongly-interacting matter. With recent high-statistics lattice QCD sim- 
ulations of the finite-temperature equation of state, a detailed comparison 
with model results becomes accessible. Such comparisons allow to estimate 
volume and truncation effects of quantities, obtained on the lattice and 
provide possible lattice extrapolation procedures to finite chemical poten- 
tial which are important to locate a critical endpoint in the QCD phase 
diagram. 

PACS numbers: 12.38.Aw, ll.10.Wx , 11.30.Rd , 12.38.Gc 

1. Introduction 

A detailed theoretical understanding of strongly-interacting matter un- 
der extreme conditions is mandatory for various heavy-ion research pro- 
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grams. The search for possible (tri)critical endpoints in the phase diagram 
is a major focus of the CBM experiment at the future FAIR facility. 

Different regimes of the QCD phase diagram can be explored by employ- 
ing various theoretical methods. Lattice QCD simulations are applicable 
at zero or imaginary chemical potentials. Most simulations for (2 + 1)- 
flavor QCD agree that chiral symmetry is restored by a smooth crossover 
transition at fi = [l] but there is still an ongoing discussion concerning 
the (pseudo) critical temperatures [2]. At finite real chemical potentials the 
fermion sign problem remains a considerable obstacle. Several extrapola- 
tion techniques to finite chemical potentials such as the reweighting method, 
imaginary chemical potential or a Taylor expansion around vanishing chemi- 
cal potentials have been proposed (for an overview see [3]). Effective models, 
such as the Polyakov-loop NJL (PNJL) or the Polyakov-loop quark-meson 
(PQM) model [3], incorporate important fundamental symmetries and the 
symmetry breaking pattern of the underlying QCD but do not suffer from 
the sign problem, large quark masses or finite volume restrictions. Further- 
more, these type of models can be used to test certain lattice extrapolation 
techniques to finite [i. For this purpose the reproduction of the lattice data 
for vanishing [/, is a basic prerequisite. 

In this talk we compare the bulk thermodynamics of several chiral 
Nf = 2 + 1 quark flavor PQM models with recent iV T = 8 lattice data 
of the HotQCD collaboration [5]. The larger quark masses used on the lat- 
tice are explicitly considered in the model comparison. Based on a novel 
differentiation technique, higher derivatives of the thermodynamic poten- 
tial can be calculated very precisely [6]. This method allows to investigate 
convergence properties of the Taylor expansion method used on the lattice. 

2. Polyakov-Quark-Meson Model 

The Polyakov-quark-meson (PQM) model for three quark flavors is based 
on the linear cr-model with quarks [7] and incorporates in addition the 
Polyakov loop field &(x). The Polyakov loop is the thermal expectation 
value of a color traced Wilson loop in the temporal direction. In the heavy- 
quark limit <I> serves as an order parameter for the confinement / deconfinement 
transition. In the deconfined high-temperature phase the center symmetry 
Z(3) of QCD is spontaneously broken and as a consequence $ is finite. 
However, with dynamical quarks the center symmetry is always broken and 
<3? as an order parameter becomes questionable. 

The PQM Lagrangian consists of a quark-meson contribution and a 
Polyakov-loop potential U(&,&), which depends on $ and its hermitian 
conjugate The uniform temporal background gauge field is coupled to 
the quarks by replacing the standard derivative in the quark contribution 
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by a covariant derivative = 
This leads to the Lagrangian 



d^i iA^ , A^ 



S^A where A^g s A"\ a /2. 



C PQ M = q(ilp-g</>5)q + Cm-U($lA],$[A}) 



(1) 



where q denotes the quark field. The interaction between the quarks and the 
meson nonets is implemented by a flavor-blind Yukawa coupling g and the 
meson matrix ^5 = Yll=o(^a/^) ( a a + il5 7T a) where the nine scalar mesons 
fields are labeled by a a and accordingly the nine pseudoscalar fields by ir a . 
The remaining, purely mesonic contribution reads 

2 / , \ 2 
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with the fields <f> = ^2 a (^a/2) (cr a + i^a)- Chiral symmetry is explicitly 
broken by the last term in Eq. ((2J) and the U (l^-symmetry by the 't Hooft 
determinant term with constant strength c. 

For the effective Polyakov loop potential Li which is constructed in terms 
of $ and <£, several implementations are available. The simplest choice is 
based on a Ginzburg-Landau ansatz [8]: 
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The cubic $ terms are required to break the symmetry of the remaining 
terms down to the center symmetry, Z(3). The potential parameters are 
adjusted to the pure gauge lattice data such that the equation of state and 
the Polyakov loop expectation values at finite temperature are reproduced. 
An improved version [2] based on the SU(3) Haar measure results in 
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a(T)$$ + b(T) In 1 - b$$ + 4 ($ 3 + $ 3 ) - 3 ($$) : 
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The parameters are again fitted to the pure gauge lattice data. The loga- 
rithmic form constrains and <1> to values smaller than one. Another choice 
invented by Fukushima [10] is 



Fuku 

rp4 
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54e~ a/T $$ + In (l - 6$$ - 3($$) 2 + 4($ 3 + $ 3 ))] (5) 



with only two parameters a and b. The parameters also result in a first- 
order transition at To ~ 270 MeV in the pure gauge sector but are not 
fitted to lattice data. This potential excludes contributions of the uncon- 
fined transverse gluons to the equation of state which are relevant at high 
temperatures [10] . This is important for the comparison to lattice data at 
higher temperatures. 
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2.1. Thermodynamic Potential 

For three quark flavors the grand potential is a function of the tem- 
perature and in general three quark chemical potentials, one for each flavor. 
Here we focus on a uniform quark chemical potential fi = fj, q = /iB/3. Since 
we consider the isospin-symmetric case with mi = m u = m^, only two order 
parameters, the non-strange a x and strange a y emerge [7]. The thermody- 
namic potential in mean-field consists of three different contributions [llj : 
the mesonic part U (o~ x , o~ y ), a quark part 0, qq and the Polyakov loop poten- 
tial 

n = U(a x ,<ry) + Slq q ((T x ,<Ty,$,$) +W ($,<&) . (6) 

The mesonic contribution has six parameters which are fitted to the vacuum. 
For example, the Yukawa coupling g is fixed to reproduce a light constituent 
quark mass of mi ~ 300 MeV. This then yields a strange constituent quark 
mass of m s ~ 433 MeV. 

The temperature- and quark chemical potential dependence of the four 
order parameters for the chiral and confinement /deconfinement transition 
are determined as solutions of the corresponding gap equations. These 
coupled equations are obtained by minimizing the grand potential, Eq. ([6|), 
with respect to the four constant mean-fields (<j) x , (°~) y i an d 

on _ on _dn _ an 

do~ ~ IkTy ~ d¥ ~ M 

where min = {a x = (cr x ) , cr y = (cr y ) ,$ = ($),$= ($)} labels the global 
minimum. 

3. Lattice comparison 

In order to compare our model results with recent HotQCD lattice find- 
ings [5] we have chosen a model parameter setup where the chiral and de- 
confinement transition temperature coincide at \x = 0. This is the case for 
m a = 600 MeV and To = 270 MeV for all used Polyakov-loop potentials [llj . 
A coincidence of both transitions around T x ~ 185 — 195 MeV is also ob- 
served in the corresponding lattice simulations [5]. However, on the lattice 
a ratio of the physical strange quark mass to the light one of 10 has been 
used which yields finally too heavy light-quark masses. For a proper com- 
parison we have therefore adjusted the pion and kaon masses in the model 
calculations accordingly and use also mx = 503 MeV and m^ = 220 MeV. 
Of course, heavier meson masses yield also slightly heavier constituent quark 
masses 

mi w 322 MeV and m s 438 MeV . 
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(a) Normalized pressure p/psB (b) Interaction measure A/T 4 

Fig. 1. The normalized pressure (left panel) and the interaction measure (right 
panel) as a function of temperature. The model calculations (PQM model with 
various Polyakov-loop potentials and the QM model) are compared to lattice data 
(N T = 8, p4 and asqtad actions) from [5]. The solid lines correspond to larger pion 
and kaon masses as used in the lattice simulations while the dashed lines denote 
the results for physical masses. 



As a consequence, higher transition temperatures, in particular for the chiral 
transition in the light sector, are found. But both transitions, the non- 
strange chiral and the deconfinement one, still coincide. The strange quark 
sector is almost unaffected 

In the model calculations all thermodynamic quantities are extracted 
from the grand potential. The pressure 

p(T, f i) = -n(T, f i) (8) 

is directly related to the thermodynamic potential with a suitable normaliza- 
tion, p(0, 0) = 0. In Fig.[T]the pressure, normalized to the Stefan-Boltzmann 
(SB) value of the PQM model, is compared to the lattice data and a quark- 
meson (QM) model calculation. As expected, the QM model [7] fails in 
describing the lattice data, while the PQM model results are in better 
agreement with the data. The best agreement is achieved with Fukushima's 
potential where the different treatment of the transverse gluons at higher 
temperatures is apparent. Around the transition the model versions with 
the polynomial or logarithmic Polyakov loop potential are closer to the lat- 
tice data. For physical meson masses (dashed lines) the pressure is smaller. 

The interaction measure, A = e — 3p, indicates the breaking of scale 
invariance and is given as a temperature derivative of the grand potential 



a A/m = _ t A (SL) 

T 4 <9T V T 4 / dT\T*J 
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In lattice simulations, this quantity can be obtained directly from the trace 
of the energy-momentum tensor 

e£ = e_-3p = A 

rpi y>4 y4 ' ^ ' 

It is more sensitive to finite volume and discretization effects than, e.g., 



the pressure (cf. [5]) which is displayed in Fig. 1(b) In the chirally broken 
phase, the lattice data are close to all model curves. Around the transition 
T ~ T x Fukushima's potential is closest to the asqtad-action data, while 
the remaining two model curves are in better agreement with the p4-action 
data. In general, the peak height decreases on larger lattices, in particular 
for the p4-action while the asqtad-action shows a weaker iV T dependence [5]. 
In contrast to the pressure, the logarithmic and polynomial Polyakov model 
version do not describe the lattice data in the symmetric phase. However, 
Fukushima's ansatz approaches the data at least up to T ~ 1.5T X . 

In addition we also show the energy density e and the entropy density 
s, which are defined as 

e = -p + Ts , s = -|p> (11) 

in Fig. [2] as a function of temperature. Similar to the pressure, Fukushima's 
potential model version comes closest to the lattice data for T < 1.5T X . 

4. Convergence of the Taylor expansion at finite fi 

At finite chemical potential the spectrum of the Dirac operator becomes 
complex making direct Monte Carlo simulations impossible. Several meth- 
ods have been developed to circumvent this problem and to access at least 
a region of small chemical potentials in the phase diagram, see e.g. [3]. 
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Fig. 3. The Taylor coefficients c@ to C22 obtained in the PQM model with the 
logarithmic Polyakov loop potential (T = 270 McV, m a = 600 MeV). 

One such approach is based on the Taylor expansion of thermodynamic 
quantities in powers of (n/T) around /i = 0. For the pressure it reads 



n=Q 



(12) 



with the coefficients 



Cn{T) 



1 d n (p(r,//)/r 4 ) 



n\ d {n/T) n 



(13) 



(i=0 



Note that only even coefficients contribute since the QCD partition function 
is CP-symmetric, i.e., Z(/jl) = Z[—\i). The coefficients have been calculated 
up to the eighth order by different lattice groups, for Nj = 2 see e.g. [12] 
and for Nf = 2 + 1 quark flavors, see e.g. [13]. However, higher orders still 
suffer from large errors. 

The application of the Taylor expansion method to the PQM model 
provides an opportunity to investigate and test the finite \x extrapolation 
since a direct evaluation of the thermodynamic potential in the model is 
possible. For example, the location of the critical endpoint (CEP) is exactly 
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known in the model calculation. This enables us to distinguish between 
divergences related to the CEP and to the breakdown of the expansion. 

However, a further difficulty emerges in the model analysis: in principle, 
the grand potential is known analytically but an analytic evaluation of the 
Taylor coefficients fails since the implicit temperature and chemical poten- 
tial dependence of the order parameters is only known numerically via the 
equations of motion Eq. ((7J). 

Standard methods for an evaluation of numerical derivatives are ham- 
pered by increasing errors, which become dominant in particular for higher 
derivative orders. In order to circumvent this caveat we have developed 
a novel numerical technique which is based on algorithmic differentiation 
(AD). With this method the evaluation of higher derivatives becomes fea- 
sible to extremely high precision. In fact, it is essentially limited only by 
machine precision. Details of this method can be found in Ref. [6]. 

With this method we could obtain the Taylor coefficients up to 22 nd 
order for the PQM model with the logarithmic Polyakov loop potential as 
shown in Fig. [3l The coefficients are very small far away from the transition 
temperature and they start to oscillate with increasing amplitude within a 
narrow temperature range around the transition temperature, i.e. 0.95T X < 
T < 1.05T X . This signals a complex /x singularity of the thermodynamic 
potential close to the real axis. In the chiral limit this singularity would be 
exactly on the real axis as a reflection of a real phase transition. 

By means of the Taylor expansion method it is further possible to study 
several thermodynamic quantities at small \x. One example is the pressure 
difference 



In Fig. |4]we show Ap for /i/T = 0.8 and fi/T = 1.0 for various expan- 
sion orders in comparison with the model result. The ratios [ijT are slightly 
below and above the critical value pi c /T c ~ 0.9 at the CEP. Although the 
transition is still a crossover at \xjT = 0.8, a divergence is observed in the 
Taylor expansion. The divergence becomes more prominent with increas- 
ing orders. However, a direct model evaluation of the pressure, labeled as 
'PQM' in the figure, shows a smooth behavior. The divergence in the coef- 
ficients is clearly a signal for the breakdown of the Taylor expansion, which 
is also related to the oscillations in the coefficients. This occurs already for 
fi/T < 1 and indicates the importance of the higher order coefficients. For 
temperatures just below the breakdown at T ~ 0.95T X the higher orders 
improve the observed agreement with the 'PQM' curve. A similar behavior 
is also seen for the fi/T = 1 case. At higher temperatures T ^> T x the coef- 
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Fig. 4. The pressure difference Ap/T 4 for different ratios of fx/T and orders of the 
Taylor expansion. The black line, labeled with 'PQM', is the model calculation. 
The vertical dashed line indicates the radius of convergence. 



ficients become small again and do not longer contribute. As a consequence 
the Taylor expansion can again reproduce the PQM result. 

For a deeper understanding of the breakdown of the Taylor expansion it 
is instructive to investigate its convergence radius. It can be obtained from 
the definition 

1/2 

r = lim T2n = lim 



02,, 



C2n+2 



(15) 



It is not yet known how well the radius r is estimated by r n for a finite 
value n. The estimated T24 is indicated by the dashed vertical line in Fig. HJ 
It coincides with the occurrence of the divergence observed in Ap at 24-th 
order. From the different r n we expect that the true convergence radius 
is smaller than r24 for the given ju/T ratios. For more details see |14j . 
A detailed study concerning the prospects of locating the CEP within the 
Taylor expansion will be given in |15j . 



5. Summary 

In this talk we have presented results for the bulk thermodynamics of 
a (2 + l)-flavor PQM-model with three different effective Polyakov loop 
potentials. The model parameters are adjusted at vanishing chemical po- 
tential to produce a coincidence of the chiral and deconfmement transition 
for all three Polyakov loop potentials. For a better comparison with the 
recent HotQCD lattice data we further tuned the pion and Kaon masses 
to these values which are used in the lattice simulations. All three PQM 
model versions yield a good reproduction of the lattice data for temperatures 
T < 1-5T X . Furthermore, the Taylor expansion of the pressure difference 
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for finite chemical potentials is also investigated. For this purpose, a novel 
algorithmic differentiation technique has been developed. The Taylor co- 
efficients up to 24 th order could be obtained in the PQM model for the 
first time. The knowledge of these higher Taylor coefficients allows for a 
systematic study of convergence properties of the expansion. 
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